Diamagnetic correction to the ^Be+ ground-state hyperfine constant* 
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We report an experimental determination of the diamagnetic correction to the ^Be+ ground state 
hyperfine constant A. We measured A = -625 008 837.371(11) Hz at a magnetic field B of 4.4609 T. 
Comparison with previous results, obtained at lower values of B (0.68 T and 0.82 T), yields the 
diamagnetic shift coeflicient k = 2.63(18) x 10^" T~^, where A{B) = AqX (l + kB^). The zero-field 
hyperfine constant Aq is determined to be —625 008 837.044(12) Hz. The y-factor ratio gi'/gj is 
determined to be 2.134 779 852 7(10) x 10~^, which is equal to the value measured at lower B to 
within experimental error. Upper limits are placed on some other corrections to the Breit-Rabi 
formula. The measured value of k agrees with theoretical estimates. 

PACS numbers: 32.10.Fn, 32.10.Fn, 32.60.+i, 32.30.Bv 
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I. INTRODUCTION 

Transition frequencies between hyperfine- Zeeman sub- 
levels of ground or metastable electronic states of atoms 
can in some cases be measured extremely accurately. One 
example is the ground-state hyperfine transition of "'^'^^Cs, 
which is currently the basis of the SI second [1]. Its fre- 
quency can be measured with a relative accuracy of bet- 
ter than 5x10^^^ [2] . Because of the accuracy with which 
the energy separations can be made in the ground states 
of Cs and other atoms, various small contributions to the 
energies can be observed and compared with calculations. 

The relative energies of the sublevels of an atom with 
electronic angular momentum J and nuclear spin I, in 
a fixed magnetic field B, are, to a good approximation, 
determined by the effective Hamiltonian 



H = hAI 3 m B fij B 

= hAI ■ J + g/fiBl ■ B + qj^j-bJ ■ B. 



(1) 



On the right-hand side of Eq. (1), the first term is the 
magnetic dipole hyperfine interaction, the second is the 

I 



nuclear Zeeman interaction, and the third is the elec- 
tronic Zeeman interaction. Here, h is the Planck con- 
stant, A is the magnetic dipole hyperfine constant, and 
fij and fij are the nuclear and electronic magnetic mo- 
ment operators, respectively. The g-factors are defined 
by gj = -Hj/iJus) and gi' = ~iii/{Ihb), where hb is 
the Bohr magneton. The prime in gi' is to distinguish 
it from the alternative definition gi = iij/{IfiN), where 
/iAT is the nuclear magneton. For I > 1 and J > 1, other 
terms, such as the electric quadrupole hyperfine interac- 
tion, should also be included on the right-hand side of 
Eq. (1). 



For J = i , the eigenvalues of H are given analytically 
by the solutions of quadratic equations. The expression 
for the energy eigenvalues is known as the Breit-Rabi for- 
mula [3, 4]. (In its original form, the Breit-Rabi formula 
did not include the nuclear Zeeman interaction, since it 
was considered to be negligible [3].) The Breit-Rabi for- 
mula for the energies of the {F, rap) sublevels in a state 
with J = i and / > i is 
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For the special case F = / + i, = ±{I + i), 
E{F,mF) = hA 
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Here, F — J + I, and rrii? is the eigenvalue of F^. At 
finite -B, the energy eigenstates are not eigcnstates of 
F^, except for the F = 1 -\- \, mp = ±(J' + 5) states. 
Nonetheless, we label them by the value of F that is 
valid at B 0. Here X = MsS(.gj - gi')/[{I + l/2)hA] 
is a dimensionless quantity proportional to B, and 7 = 
gi'/gj is the g-factor ratio. The ± sign in Eq. (2) cor- 
responds to the states labeled by = / ± i. For the 
ls^2s^Si/2 ground electronic state of ^Be+, which was 
the subject of this study, ^ = | . The value of gj for the 
ground electronic state of ^Be+ has been determined by 
measuring the ^Be"*" cyclotron frequency and a hyperfine- 
Zeeman transition frequency at the same magnetic field 
[5]. The value is gj = 2.002 262 39(31), calculated with 
the use of the best current value of the proton-electron 
mass ratio [6]. For ^Be+, X « -22.414B(T). 

There arc several ways in which the experimental en- 
ergy separations can deviate from those predicted by the 
Breit-Rabi formula. For / > i, it is possible that, at a 
fixed value of no values of the parameters A, g/ /gj, X 
can be found that are consistent with all of the measured 
energy separations. It is also possible that the values 
of the parameters determined at one value of B are not 
consistent with those determined at another value. 

There arc several possible sources of deviations from 
the Breit-Rabi formula. If there is another electronic 
energy level that is close in energy, hyperfine or Zeeman 
interactions can mix the electronic states. For the ground 
state of an alkali atom or an alkali-like ion, there are no 
nearby electronic states, so such effects are small. More 
important are diamagnetic contributions to the interac- 
tion between the atom and the magnetic field that are 
neglected in the effective Hamiltonian given by Eq. (1). 
Diamagnetic corrections to the Breit-Rabi formula were 
first considered by Bender [7]. He calculated the size of 
the deviation in the ground state of ""^^^Cs to be equiva- 
lent to a fractional shift in A of SA/A = 3.9 x 10"^° B^, 
where B is expressed in teslas. This effect, called the 
dipole diamagnetic shift in atomic hyperfine structure, 
is due mainly to a magnctic-field-induced change in the 
electronic spin density at the nucleus. 

Measurements of magnetic-field-dependent deviations 
from the Breit-Rabi formula in the ground state of Rb 
were made by Larson and coworkers [8 10]. The dipole 
diamagnetic shift was observed experimentally in the hy- 
perfine structure of ^^Rb [8] and later in ^'^Rb [9]. A 
quadrupole diamagnetic shift was observed in ^'"'Rb and 
^^Rb [9]. In contrast to the dipole shift, the quadrupole 
shift can be thought of as a magnetically induced elec- 
tric quadrupole hyperfine interaction, which would be 
absent in a pure J = ^ state. The diamagnetic poten- 
tial, which contains a rank-2 spherical tensor part, breaks 
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the spherical symmetry, so that the electronic state is no 
longer an exact eigenvalue of J^. The signature of the 
quadrupole diamagnetic term is an energy shift propor- 
tional to [/(/ +1)- 3mj]QB'^/[I{2I - 1)], where Q is 
the nuclear quadrupole moment. In Rb, the quadrupole 
shift is smaller than the dipole shift by about three orders 
of magnitude. Another magnetic-field-dependent energy 
term was observed in ^^Rb and ^''Rb [10]. The term 
was explained by Fortson [11] and is called the hyperfine- 
assisted Zeeman shift [12]. The shift of a level is propor- 
tional to [m/^mj — /(/ -|- l)mj -|- mi/2]{gj')'^B and is 
due to mixing of higher electronic states with reversed 
electronic spin into the ground electronic state by the 
magnetic dipole hyperfine interaction. 

The ground-state hyperfine constant of ^Be+, A, was 
measured with a fractional uncertainty of 2.4 x 10~^ by 
Vetter et al. by rf-optical double resonance [13]. The 
fractional uncertainty of A was decreased to 1.6 x 10~^^ 
by Wineland et al., in measurements made with laser- 
cooled ions in a Penning trap [5]. The low uncertainty 
was due mainly to the use of transitions for which the 
first derivative of the frequency with respect to B is zero. 
Nakamura et al. measured A with a fractional uncer- 
tainty of 1.2 X 10~^ in laser-cooled ions in a linear rf 
trap, at S = 0.47 T [14]. Their value of A differed from 
that of Ref. [5] by about two standard deviations. Okada 
et al. [15] have measured A for ^Be''" in a linear rf trap. 

Based on theoretical considerations and the experi- 
mental results for Rb, the ground-state hyperfine con- 
stant of ^Be+, A, is assumed to have a weak quadratic 
dependence on B such that A{B) ^ Aq x {1 + kB^). 
Transition frequencies measured at different values of B 
are used to determine the diamagnetic shift coefficient k. 

The present experiment, on the measurement of the 
hyperfine- Zeeman transition frequencies in the ground 
electronic state of ^Be"*" in a high magnetic field {B 
= 4.4609 T), is described in Sec. II. In Sec. HI, the 
high- field results are combined with the previous, lower- 
magnetic-field measurements to obtain a value for k. A 
theoretical estimate of k is given in Sec. IV. 



II. HIGH-FIELD EXPERIMENT 

A. Atomic energy levels and transitions 

There are three unknown variables (^4, gi' /g.j, and X) 
in Eq. (2), and the measurement of three transition fre- 
quencies in the ground state at a fixed value of B will 
determine these three variables. We experimentally de- 
termined the value of A and gi' /gj at B ps 4.4609 T by 
measuring the three transition frequencies labeled /e, /i, 
and /2 in Fig. 1. While three frequencies are enough to 
determine A and gi' /g.j, we also measured a fourth fre- 
quency /a, to check for consistency. The typical period 
required for a complete set of frequency measurements 
needed to determine A, gi'/gj, and X was 30 to 40 min- 
utes. 
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FIG. 1. Energy level structure of ''Be+ at B « 4.4609 T. 
/e « 124 GHz, /i f» 340 MHz, /a f» 288 MHz, and /g ^ 287 
MHz. The frequency tuning of the repump laser is shown for 
fast repumping of the electron spin-flip transition. For the 
nuclear spin-flip measurements the repump laser was tuned 
approximately 500 MHz lower than the cooling transition. 



We trapped fewer than or approximately 10'^ ions in 
a Penning-Malniberg trap and cooled them to approxi- 
mately 1 mK by Doppler laser cooling. The cooling laser 
also optically pumped the ions into the {rrij, mj)=(|, i) 
state, labeled as the initial state \i) in Fig. 1. [Here the 
states are labeled by the (rn/, mj) quantum numbers of 
their largest components.] The basic experimental pro- 
cedure for measuring the different transition frequencies 
was to (1) turn off the cooling laser, (2) probe the desired 
transition with the appropriate rf or microwave radiation, 
(3) measure the population of the ions remaining in \i) 
with the fluorescence induced by the cooling laser, (4) 
repump all ions to \i) with the cooling laser and an addi- 
tional repumping laser. We used the same ions repet- 
itively to measure all transition frequencies. We first 
discuss the basic experimental setup and the 124 GHz 
microwave system. We then discuss in more detail the 
measurements of the different transition frequencies and 
the determination of A and gi' / gj at high magnetic field. 



B. Experimental setup 

1. Penning trap 

Figure 2 shows a sketch of the Penning trap used for 
the high-_B measurements. The trap and the basic ex- 
perimental setup have been described previously [16- 
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FIG. 2. Schematic diagram of setup. Figure is not to scale. 
The trap diameter is 4 cm. The electrodes used to apply the 
rotating wall field are not shown. The direction of the side- 
view light collection and the direction of the perpendicular 
cooling beam form a 60° angle in a plane perpendicular to 
the magnetic field axis. A side-view image of a plasma with 
approximately 500 ions is shown. The diameter of the fiu- 
orescing Be^ ion plasma is 400 /^m. Heavier-mass impurity 
ions are located at larger radii than the ®Be^ ions [16]. 



18]. The 4.4609 T magnetic field of a superconducting 
solenoid with a 125 mm room-temperature bore produces 
a ^Be"*" cyclotron frequency of Vic — 2tt x 7.602 MHz. 
The long-term drift of the magnetic field was less than 
one part in 10^ per hour, resulting in an average drift 
of /c of less than 3 kHz per hour. Magnetic field shifts 
due to changes in the magnetic environment (due, for 
example, to movement of Dewars or to other activity 
in neighboring laboratories) were reduced by collecting 
data during the night. We did not stabilize the pressure 
in the magnet Dewar or actively cancel external mag- 
netic field noise, which likely would have improved the 
long-term magnetic field stability [19]. The field of the 
superconducting magnet was found to have fiuctuations 
which were fast compared to 20 Hz, superimposed on 
the slow magnetic field drift and noise. The frequency 
spectrum of the fast fluctuations contained a continu- 
ous part and narrow peaks between 30 and 300 Hz [18]. 
The integrated noise of the fast fluctuations produced 
5B/B K. 10~^ variation in the magnetic field for mea- 
surements separated by greater than 0.1 s. The fast fiuc- 
tuations contributed to the linewidth and coherence of 
the electron spin- flip measurement, but had no signifi- 
cant impact on the nuclear spin- flip measurements (/i, 
/2, and /s). (The fast fluctuations only produce a phase 
modulation of a few milliradians on the nuclear spin-flip 
transitions.) Recent work, which will be discussed in 
a separate publication, indicates that the fast fluctua- 
tions are fluctuations in the homogeneous field produced 
by the superconducting magnet, which can be mitigated 
through vibration isolation of the magnet. 

The Penning trap electrode structure consists of a 
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stack of four cylindrical electrodes. The inner diame- 
ters of the cylinders are 4.1 cm and the combined length 
of the four cylinders is 12.7 cm. We typically oper- 
ated the trap with the central cylindrical electrodes (the 
"ring" electrodes) biased at —500 V and the outer cylin- 
drical electrodes (the "endcap" electrodes) grounded, 
which resulted in ^Be+ single-particle axial and mag- 
netron frequencies of, respectively, = 27r x 565 kHz 
and ojm — 2ttx 21.1 kHz. 

Due to the crossed electric and magnetic fields in a 
Penning trap, an ion plasma undergoes a rotation about 
the magnetic field axis. In thermal equilibrium this rota- 
tion is rigid [20] , and we use ujr to denote the plasma ro- 
tation frequency. The rotation frequency of the ^Be+ 
plasma was precisely controlled with a rotating electric 
field (a rotating wall) [17, 21]. A rotation frequency LOr of 
2iTX 30 kHz or less was used, which produced planar plas- 
mas (oblate spheroids) like that shown in Fig. 2 with ion 
densities of approximately 8 x 10^ cm~^. The measure- 
ments presented here were obtained on small ion plasmas 
of fewer than 10^ ^Be"*" ions. The small axial extent of the 
plasmas (typically less than 50 /um) reduced the effect of 
axial gradients in the magnetic field. Axial magnetic field 
gradients were shimmed to be less than two parts in 10* 
per mm, which resulted in an axial magnetic field inho- 
mogeneity of less than one part in 10^ over a 50 /im axial 
extent. We found no evidence for any inhomogeneous 
broadening of the different resonance curves discussed in 
Sees. HC and HD. 



2. Laser cooling, state preparation, and detection 

Doppler laser-cooling was carried out on the 313 nm 
25'i/2 (to/ = |,TOj = i) -> 2p3/2 (rni = |,mj = |) 
transition (see Fig. 1). The 313 nm light was gener- 
ated by frequency-doubling the output of a dye laser 
at 626 nm. The axial and perpendicular cooling beams 
cooled the motion parallel and perpendicular, respec- 
tively, to the magnetic-field axis [16, 22]. The axial cool- 
ing beam had a 1 mm waist diameter, a power of approxi- 
mately 1 mW, and a polarization that was either linear or 
circular (o"^). The axial cooling beam was aligned with 
the magnetic field axis to better than 0.01°. The perpen- 
dicular cooling beam was linearly polarized in a direction 
perpendicular to B, focused to a waist diameter of ap- 
proximately 50 /im, and had a power of approximately 1 
/uW. A double-pass acousto-optic modulator was used to 
rapidly switch the cooling beams off (in less than 1 /xs) 
before applying rf or microwave radiation to drive the de- 
sired ground-state transitions. The cooling beams were 
switched back on after the rf or microwave radiation was 
switched off. 

The population of the "^81/2 [mi — ^,mj — ^) state 
(the initial state \i) in Fig. 1) was measured through 
the cooling-laser-induced resonance fluorescence. An //5 
imaging system was used to image the ^Be"*" ion fluores- 
cence onto the photocathode of a photon-counting imag- 



ing tube (quantum efficiency was approximately 5 %). 
The total imaging tube count rate was proportional to 
the \i} state population. The total photon count rate 
was recorded for 0.5 s both before and after applying the 
rf or microwave radiation. The ratio of these two count 
rates (with small corrections for repumping effects) mea- 
sured the fraction of the ions remaining in the \i) state. 

The cooling radiation optically pumped more than 
94 % of the ions into the ^51/2 (f , 5) state, i.e., the lower 
level of the cooling transition [23, 24]. This was a non- 
resonant optically pumping process with a time constant 
of approximately 5 s for the cooling laser parameters in 
this experiment. The repumping time on the electron 
spin- flip transition (/e) was reduced to less than 1 ms by a 
second frequency-doubled dye laser (labeled the repump 
laser in Fig. 1). The repump laser was turned on after the 
second 0.5 s detection period (the detection period after 
the applied rf or microwave radiation). The /i, /2, and 
fs transitions involve a change in the ^Be+ nuclear spin 
orientation. For example, in /2 the nuclear spin changes 
from TO/ = I to TO,/ = i. Optical repumping back to ji) 

occurred through the ^5*1/2 (5, ^) — > ^^'3/2 (^i |) transi- 
tion and the small admixture of different (to/, toj) states 
in the 2p'^P^i2 manifold [24]. We reduced this repump- 
ing time somewhat by tuning the repump laser frequency 
between 400 MHz and 600 MHz below that of the cooling 
transition. This had the added benefit of maintaining a 
cold-ion plasma when most of the ions arc driven to the 
(to/ = ^,TOj = i) state. Presumably this was because 
the frequency of the repump laser was now below that 
of the "^81/2 (|)|) transition. A simi- 

lar improvement in the repumping and plasma stability 
was also observed for the /i and /s transitions with a 
"far-detuned" laser tuned 400 MHz to 600 MHz below 
the cooling transition. The repump or far-detuned beam 
was directed along the magnetic fleld axis of the trap, as 
shown in Fig. 2. The beam waist diameter was approxi- 
mately 0.5 mm; the power was a few milliwatts. 



3. Microwave apparatus 

A sketch of the 124 GHz microwave system used to 
measure /e is shown in Fig. 3. Reference [18] further 
discusses the microwave system and its use in quantum 
information experiments. A Gunn diode oscillator gener- 
ated 30 mW of microwave power, and its frequency was 
coarsely set to approximately 124 GHz by a manually 
tuned microwave cavity. The microwave radiation was 
transmitted through WR-8 wave guides and launched to 
free space through a pyramidal rectangular microwave 
horn. A small fraction of the microwave power (—10 dB) 
was mixed with the 8th harmonic of a 15.5 GHz dielec- 
tric resonator oscillator (DRO). The intermediate fre- 
quency (IF) signal from the harmonic mixer was sent to 
a phase-locked loop (PLL) controller and phase-locked to 
a 76 MHz reference frequency generated by direct digi- 
tal synthesis (DDS). The microwave frequency and phase 
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were controlled by changing the frequency and phase of 
the DDS signal. The DDS was controlled by computer 
through a parallel interface. All of the frequency syn- 
thesizers used in the experiment were referenced to the 
same passive hydrogen maser, including the DRO and 
the DDS. The frequency of the passive hydrogen maser 
was calibrated relative to that of the NIST atomic time 
scale. 



spectrum 



analyzer 



oscilloscope 



100 kHz 



microwave 
horn 

switch 
c TTL 



76 MHz 



harmonic mixer 







IF mon IF in 

Phase-Lock 



Loop 



ref in 



Gunn 
out 



IF 



RF 
L0x8 



15.5 GHz 



profile 
switch 
TTL^ 



\d-20 dB 
i76 MHz 



DDS 



DRO 



100 MHz 



1 GHz 



cioctc 
generator 



10 MHz 







(Z) 



attenuator 



1 dB coupler 



phase shifter 

124 GHz 
Gunn diode 



H maser 



FIG. 3. Schematic diagram of the 124 GHz microwave source. 
A Gunn diode oscillator generates 30 mW of microwave power 
that was transmitted through WR-8 wave guide (shown in 
gray) and emitted to free space via a pyramidal microwave 
horn. The frequency of the Gunn diode oscillator was phase- 
locked to a signal derived from a hydrogen-maser frequency 
standard. 



The microwave radiation emitted from the horn was 
rapidly switched on and off with a reflective PIN diode 
switch. The ratio of the high and low power states was 
26.5±0.9 dB. Switching between the high- and low-power 
states could cause the PLL to lose phase-lock because of 
the change in the reflected signal. To avoid losing phase- 
lock, we added a phase shifter between the Gunn diode 
oscillator and the PIN diode switch. (The large fringing 
field of the magnet made use of an isolator impractical.) 
The phase shifter required careful adjustment to achieve 
a condition where the Gunn diode oscillator would not 
lose lock when the microwave power was switched. 

We used quasi-optical techniques (with a horn and 
Teflon^ lens) to couple the microwave radiation to the 
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FIG. 4. Schematic diagram of the quasi-optical coupling of 
124 GHz microwave radiation to the ions. Figures are not to 
scale, (a) Side view, (b) Top view. The inner diameter of the 
super-conducting magnet was 12.5 cm, the distance from the 
horn to the lens was 29 cm, and from the lens to the ions is 
28 cm. The horn was shifted off-axis in order to avoid hitting 
the horn with the axial laser beams. The lens was also shifted 
about half this amount to focus on the ions. The lens diameter 
was 10.2 cm. It had a cut on the side to make a room for the 
//5 optics used for the side-view camera. Electron spin- flip 
TT-pulse periods of 100 /xs were obtained with this setup. The 
position of the rf antenna is also shown. 



ions, as schematically shown in Fig. 4. The pyramidal 
horn coupled mainly to the Gaussian TEMOO mode with 
an initial waist diameter of approximately 0.7 cm. The 
waist diameter w is defined by P{r)/P{0) — e^^^"^/^^ , 
where P{r) is the power per unit area at a radius r. The 
waist diameter of the Gaussian beam increases as it prop- 
agates. The beam was focused to a waist of about 0.7 
cm at the ions by the hyperbolic surfaces of the lens [25] , 
which was 29 cm below the horn and 28 cm above the 
ions. The lens had a cut near the side for the //5 side- 
view optics and a 4.8 mm diameter hole near center to 
pass laser beams (see Fig. 4). The horn was shifted off 
axis to avoid the laser beams, and the lens was shifted 
accordingly to center the microwave focus on the ions. 
The microwave system was on an x-y mechanical stage, 
and the position of the horn was adjusted to maximize 
the coupling of the microwave radiation with the ions. 

Electron spin-flip 7r-pulse periods as short as 100 
were obtained with this microwave system. 7r-pulse fideli- 
ties of better than 99.9 % were measured with random 
benchmarking on plasmas consisting of a single plane 
[18]. The Ramsey free-induction decay (that is, the free- 
precession period in a Ramsey experiment where the 
fringe contrast has decayed by 1/e) was measured to be 



^ Teflon is a registered trademarli of ttie Dupont Company. Men- 
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that this material is endorsed by NIST or that it is necessarily 
the best material for the purpose. 
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T2 « 2.4 ms, limited by the fast magnetic field fluctua- 
tions. 



C. Electron spin-flip measurement 

Figure 5 shows an electron spin-flip resonance obtained 
with a 600 /iS square Rabi pulse. The data were fitted to 
the expected Rabi resonance curve [4], 



(26)2 



(/e - fr + m- 



■ sm 



(4) 

Here Pi is the probability of an ion to be in state h = 
f2/27r, where Q, is the Rabi frequency, t — 600 /is is the 
microwave pulse duration, / is the microwave frequency, 
and /e is the electron spin-flip resonance frequency. From 
the fit to the data in Fig. 5, we determine a value for the 
electron spin-fiip frequency /c = 124 076 860 036 ±15 Hz. 
The uncertainty obtained from the fit we define to be the 
internal error, and for electron spin- flip resonance curves 
taken under conditions similar to that shown in Fig. 5, 
the internal error was typically less than 20 Hz. Because 
/e is roughly proportional to B, 20 Hz corresponds to a 
1.6 X 10^^° fractional measurement of B, a reduction by 
about a factor of five, due to averaging, from the shot- 
to-shot variation in B. 
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FIG. 5. Electron spin-flip Rabi resonance. Tiie microwave 
power is adjusted to produce a 7r-pulse on resonance in 600 //s. 
The solid line is a fit to the expected Rabi resonance, as 
discussed in the text. All experimental measurements were 
equally weighted in the fit. 



The resonance curve in Fig. 5 took a few minutes to 
obtain. If we took many resonance curves over a longer 
period, the scatter in the fitted values for /e was larger 
than the internal error of an individual fit, due to slow 
drift and fluctuations in B. We define the external error 



to be the standard deviation of /c, determined from many 
separate scans. The typical external error for data taken 
over a 30-minute period was greater than 40 Hz. We 
found the external error to be smallest between 10 PM 
and midnight local time. Since the measurement of /o 
was limited by the stability of B, reducing the internal 
error by narrowing the line width with longer Rabi pulses 
would not have benefited us. 

To determine the hyperfine constant A, we cycled be- 
tween measurements of /o and measurements of the nu- 
clear spin-fiip frequencies discussed in the next section. 
To minimize the effect of drift and slow fluctuations in B, 
it was important to complete one cycle of measurements 
as rapidly as conveniently possible. We found that we 
could complete a cycle of measurements more rapidly if 
we did not change the frequency of the repump laser to 
the repumping transition (see Fig. 1) for the /e measure- 
ments. Therefore, we set the frequency of the repump 
laser to the far-detuned position (400 MHz to 600 MHz 
lower than the cooling transition) for all measurements 
of /c, /i, /2, and /s. 



D. Nuclear spin-flip measurements 

The rf radiation used to drive the nuclear spin-flip 
transitions was generated by mixing the output of an 
80 MHz synthesizer having 1 mHz resolution with a 
higher-frequency synthesizer that had lower frequency 
resolution. For measurements of J2 and /a (approxi- 
mately 290 MHz), the higher- frequency synthesizer was 
set to 220 MHz. For measurements of /i (approximately 
340 MHz) the higher-frequency synthesizer was set to 
280 MHz. Switching of the rf was done with a switch 
having approximately 90 dB isolation. The rf radiation 
was coupled to the ions through a two-turn rf loop an- 
tenna, placed near ions, outside the vacuum envelope (see 
Fig. 4). 

Figure 6 shows an /2 resonance curve obtained with the 
Ramsey method. The rf power was adjusted to achieve a 
7r/2 pulse in 0.5 sec. The two 7r/2 pulses were separated 
by 4 s. After the Ramsey sequence, both the cooling laser 
and the far-detuned laser were turned on simultaneously. 
The power of the far-detuned beam was adjusted so that 
the dark state repumped with a 1/e time constant of 
approximately 2 s. To avoid any significant ac Zeeman 
shifts due to the finite isolation (26 dB) of the microwave 
switch, the microwave frequency was detuned from res- 
onance with /c by 1 MHz during the /2 measurement. 
Fitting the data of Fig. 6 to a sinusoidal curve gives /2 — 
288 172 932.435 3(7) Hz. The typical external error from 
measurements taken over a 30 to 40 minute measurement 
cycle was approximately 5 mHz. A 5 mHz external error 
with the 6.5 kHz/mT sensitivity of this transition implies 
a fractional magnetic field stability of 2 x lO"""^" over a 
typical 30 to 40 minute period, which is comparable to 
what was observed on the electron spin-flip transition. 

Figure 7 shows an /i resonance curve obtained with 



7 



Doppler cooling laser 



0.5 s 0.5 s 




0.0 0.1 0.2 0.3 0.4 



rf frequency - 288 172 392.220 (Hz) 

FIG. 6. /2 resonance obtained with a 4 s Ramsey free- 
precession period. The rf power was adjusted to apply a 
7r/2-pulse in 0.5 s. 



Doppler cooling laser 




0.0 0.1 0.2 0.3 0.4 



rf frequency - 339 961 281 .72 (Hz) 

FIG. 7. An /i resonance obtained with a 4 s Ramsey free- 
precession period. The rf power was adjusted to apply a 
7r/2-puIse in 0.5 s. The microwave frequency was shifted by 
±100 kHz during Ramsey sequence. 



the Ramsey method. We first transferred the popula- 
tion in |i) to the (m/ — ^,mj — — i) state with a 600 
/iS TT-pulse. This was followed by the Ramsey interroga- 
tion, as shown in Fig. 7. We then used a second 600 fis 
microwave 7r-pulse to transfer any ions remaining in the 
(m/ = fjTOj = — ^) state to The ion population in 
\i) was then detected by the laser-induced fluorescence. 



We shifted the microwave frequency by ±100 kHz from 
resonance during the Ramsey interrogation of /i. This 
prevented driving the /c transition with microwave ra- 
diation that leaked through the microwave switch. The 
DDS frequency could be switched by as much as 100 kHz, 
and we could still keep the Gunn diode oscillator phase- 
locked. A 100 kHz offset produced a 2 mHz ac Zeeman 
shift due to the microwave leakage through the switch. 
A measurement of /i consisted of taking two scans with 
alternate signs of the microwave detuning. The /i tran- 
sition frequency was determined by fitting the average 
of the two scans. A fit to the data in Fig. 7 provides 
/i = 339 961 281.917 0(8) Hz. The external error from 
measurements taken over a 30 to 40 minute period was 
typically 5 mHz, about the same as for the /2 measure- 
ments. 




0.0 Lj . 1 . 1 . 1 . J 

0.0 0.1 0.2 0.3 0.4 

rf frequency - 286 586 653.95 (Hz) 



FIG. 8. An /s resonance obtained with the Ramsey method. 
The power of the rf was adjusted to apply a 7r/2-pulse in 
0.5 second. 



Figure 8 shows an /a resonance curve obtained with a 
4 s Ramsey free-precession period. The pulse sequence 
was very similar to the one for /i. The microwave tt- 
pulse in the /i measurement was replaced with a 0.5 s 
rf TT-pulse to transfer the population from \i) to the 
(to/ = ^,TOj = i) state. Then the rf frequency and 
amplitude were changed, and we applied a Ramsey se- 
quence with a 4 s free-precession period. The frequency 
of the microwave radiation was detuned from resonance 
by 1 MHz during the /s measurement. A fit to the data 
in Fig. 8 gives /g = 286 586 654.158 5(14) Hz. 

The length of the Ramsey free-precession periods in 
the nuclear spin-flip measurements (4 s) was limited by 
the heating, presumably due to collisions with residual 
gas molecules, that occurred when the cooling laser was 
turned off [16]. For free precession periods longer than 
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4 s, the ion fluorescence decreased, due to the increase in 
the Doppler width of the coohng transition. This added 
noise and complicated the signal analysis. Longer free- 
precession periods (20 s to 100 s) have been used with 
sympathetic cooling in previous low-_B measurements [26, 
27]. The nuclear spin-flip measurements at 4.4609 T were 
limited by magnetic-field instabilities, so there was no 
compelling reason to use longer free-precession periods. 



E. Experimental results 

To determine the hyperfine constant A, resonance 
curves such as those shown in Figs. 5-7 were taken in 
succession, as shown in Fig. 9, and fits to the resonance 
curves were used to determine /o, /i, and f2- A typi- 
cal measurement cycle consisted of an /o measurement, 
followed by an /i measurement, followed by another /c 
measurement, followed by an /2 measurement, followed 
by a final /e measurement. One measurement cycle took 
30 to 40 minutes to complete. The average of the three 
/e measurements in one cycle was use to determine /c. 
The uncertainty in /e was taken to be the external error 
from the scatter in the three measurements, which was 
typically about 40 Hz. The uncertainties assigned to /i 
and /2 were the external errors from the scatter in the 
/i and /2 measurements from consecutive measurement 
cycles, which was typically about 5 mHz. 

All known systematic errors in the nuclear spin-flip res- 
onance frequency measurements, other than those due to 
the magnetic field instability, were less than 1 mHz. The 
largest systematic error is due to microwave radiation 
leaking through the microwave switch. As discussed in 
Sec. H D, this produced a 2 mHz shift in the /i resonance 
curve. However, by taking data with the microwave 
frequency shifted off resonance by both -1-100 kHz and 
— 100 kHz, this shift could effectively be canceled. During 
a 4 s nuclear spin-flip measurement, the ion temperature 
increased due to collisions with the room temperature 
residual background gas. Previous studies indicated that 
the temperature increase over a 4 s period is limited to 
a few kelvins [16, 22]. However even a 10 K temper- 
ature would produce only an approximately —0.1 mHz 
time-dilation shift in the measured ^Be''" nuclear spin- 
flip frequency. We performed some simple checks for un- 
known systematic errors by varying the length of the Rabi 
pulse in the fe measurement and the length of the free- 
precession period in the nuclear spin-flip measurements. 
In addition, we took some /i and /2 measurements with 
an 8 s Rabi pulse. No systematic dependencies were ob- 
served at the level permitted by the magnetic field sta- 
bility. 

For each measurement cycle, the Brcit-Rabi formula 
[Eq. (2)] was used to solve for values of A, gi' /gj, and 
X. The uncertainties in these values were determined 
by using the Breit-Rabi formula to solve again for A, 
gi'/gj, and X, but with /e, /i, and /2 set to the lim- 
its of their uncertainties. We conservatively assigned the 




23:00 00:00 01:00 02:00 

local time (hour:minute) 



FIG. 9. (Color online) Summary of measurements taken 
during a 4 hour run at night. The measured electron 
and nuclear spin-flip frequencies (relative to fixed frequen- 
cies /oo = 124 076 861270 Hz, /lo = 339 961281.922 Hz, 
/20 = 288 172 932.440 Hz, and /30 = 286 586 654.160 Hz) are 
plotted against the time the measurement was made. The 
measurements are grouped into cycles, labeled 1-6 on the 
graph. Each cycle was used to determine a value of A and 
gi' /gj- The error bars are the internal errors obtained from 
the fits to the resonance curves. Due to magnetic field drift, 
/e drifted down by about 100 Hz during this run. This is 
consistent with the drifts observed with /i and /2, where /i 
decreases and /2 increases with decreasing magnetic field. 



largest uncertainty that could be obtained from the dif- 
ferent combinations of limits. For example, if Sfc, Sfi, 
and 6/2 arc the uncertainties in /c, /i, and /2, solving 
the Breit-Rabi formula with the frequency values fo+Sfc, 
fi + <5/i, and /2 -I- 6/2 gives the largest uncertainty for A, 
while solving the Breit-Rabi formula with the frequency 
values fc + Sfc, fi - Sfi, and /2 + (5/2 gives the largest 
uncertainty for gi' /gj- 

Figure 10 summarizes the measurements of A at high 
magnetic field. Four different sets of data were taken on 
four different dates over a period of two months. Each set 
of data consisted of at least two and as many as six mea- 
surement cycles. The consistency of the data is good. 
The standard deviation from the scatter of the 15 dif- 
ferent measurements of A is 7 niHz, which is slightly 
less than the average 11 mHz uncertainty of an individ- 
ual measurement. For the determination of A we use a 
weighted average of the data shown in Fig. 10 and con- 
servatively assign an 11 mHz uncertainty, the average 
uncertainty for a single measurement cycle. An 11 mHz 
uncertainty corresponds to about a 3 x 10^^" fractional 
magnetic field instability, about a factor of three below 
the shot-to-shot fluctuations in the magnetic field. We 
believe that an assignment of a smaller uncertainty would 
require a careful study of the statistics of the magnetic 
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FIG. 10. Summary of the measurements of the high-field 
hyperfine constant A. Each of the four runs is a series of 
determinations of A made on the same date. Each A value is 
determined from the Breit-Rabi formula and one cycle of /e, 
/i , and /2 measurements. The assignment of uncertainties is 
described in the text. 



The energy shift of an {mj,mj) state due to the 
quadrupole diamagnetic shift, at a fixed value of B, has 
the form [9] 



Eq = hfgB 



3to|] 



/(27-1) 



(7) 



If the only correction to the Breit-Rabi formula is given 

by Eq. (7), then agreement of the measured and predicted 
values of f-^ to less than 10 mHz sets a limit |/q|-B^ < 
5 mHz at 's = 4.4609 T, or |/q| < 2.5 x lO""* Hz T'^. 

The hyperfine-assisted Zeeman shift of an (m/,TOj) 
state has the form [11] 

Ehz = 2hpiizB[mjmj - 1(1 + l)mj + m//2]. (8) 

If the only correction to the Breit-Rabi formula is given 
by Eq. (8), then agreement of the measured and predicted 
values of /a to ±10 mHz sets a limit |2/3hz-B| < 6.7 mHz 
atB = 4.4609 T, or |/3hz| < 7.5 x lO""^ Hz T-\ 



III. DETERMINATION OF fe FROM 
HIGH-FIELD MEASUREMENTS AND 
PREVIOUS LOW-FIELD MEASUREMENTS 



field fluctuations. The result is 

A(4.4609 T) = -625 008 837.371(11) Hz. (5) 

Although A was the main focus of this study, a value of 
the fi(-factor ratio is also obtained from the same analysis: 

gj'/gj = 2.134 779 852 7(10) x 10"^. (6) 

In four of the data cycles shown in Fig. 10, the fourth 

frequency /s was measured in addition to /o, /i, and /2. 
These measurements were used to place upper limits on 
corrections to the Breit-Rabi formula at a fixed value of 
B. The results are listed in Table I. If the Breit-Rabi 
formula is assumed to be correct, and A and g/ /gj are 
fixed at the values given by the complete set of 15 data 
cycles [Eqs. (5) and (6)], then any of the three frequen- 
cies /e, /i, or /2 can be be used to determine X for a 
given data cycle. This value of X can then be used to 
predict /a for that data cycle. In practice /2 tended to 
yield the most consistent values of X. The rms differ- 
ence between the measured value of /a and the value 
predicted from the measurement of /2 was 7 mHz. This 
is consistent with the noise expected from magnetic field 
fluctuations. We use these results to place an upper limit 
of 10 mHz on any shifts of /a at B = 4.4609 T due to 
corrections to the Breit-Rabi formula. More specifically, 
we can set limits on corrections to the Breit-Rabi formula 
having the form of the quadrupole diamagnetic shift or 
of the hyperfine-assisted Zeeman shift. Since all of the 
measurements were made at nearly the same value of B, 
this test is not sensitive to modifications to the Breit- 
Rabi formula that amount to a dependence of either A 
or gi'/gj on B. 



The present experimental rcstilts can be combined with 
previous measurements made by some of the present au- 
thors at lower values of B to determine the S-dependence 
of A or gi' / gj- Preliminary values of A and gi jgj were 
given in Ref. [5] but not the transition frequencies on 
which they were based. We now supplement Ref. [5] with 
the transition frequencies and a final determination of A 
and gi jgj- Two nuclear spin-flip frequencies, labeled 
1 and 3 in Table II, were measured near two particular 
values of B where the first derivatives of the frequencies 
arc zero. Electron spin-flip frequencies, labeled 2 and 4 
in Table II, were measured at the same two values of B. 
The experimental methods have been described in detail 
[5, 26-28]. Transition 3 is known better than transition 
1 because it was studied for use as a frequency standard 
[27, 28]. The value we report for transition 3 in Table II 
is slightly different than that reported in [28] because it 
includes additional measurements made in 1988 and 1989 
as well as an evaluation of the background pressure shift. 

If no account is taken of any B-dependence of A or 
gi' I g.h the four frequency measurements, together with 
the Breit-Rabi formula, yield a system of four equations 
with four unknowns, A, gi jgj^ X\^ and X2, where X\ 
is the value of X for transitions 1 and 2, and Xi is 
the value of X for transitions 3 and 4. The frequencies 
given in Table II yield A = -625 008 837.053(11) Hz and 
gj' jgj = 2.134 779 8518(23) x lO"*. The precise values 
of X\ and Xi are not important, since they reflect only 
the value of B at which the experiment was performed, 
not any intrinsic property of the ^Be"*" ion. Comparing 
these results to Eqs. (5) and (6), we see that there is 
clear evidence for i3-dependencc of A, but that gi jgj is 
independent of B to within experimental error. 
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TABLE I. Measured values of /a and values of /a predicted from the Breit-Rabi formula. All frequencies are in hertz. 

/2 (measured) fs (predicted from /2) /a (measured) fs (predicted) — /a (measured) 

288172 931.6510 286 586 653.404 9 286 586 653.400 3 00046 

288172 931.609 7 286 586 653.365 8 286 586 653.372 7 -0.0069 

288172 932.447 286 586 654.156 8 286 586 654.1671 -0.0103 

288172 932.453 5 286 586 654.162 9 286 586 654.157 8 0.0051 





TABLE 11. Low-B 


resonance frequencies used to determine A and gi' /gj- 




Label 


{mi,mj) -n- (mi ,mj') 


B(T) 


Frequency (Hz) 


Uncertainty (Hz) 


1 

2 
3 
4 


( 2^ ~ 2 ) ^ (2^ ~ 2 ) 

(2^' ^^2) ^ (2^ 2) 
("^2' 2) *^ ^~^3 ^1 
(~2'~2) ^ (~2'2) 


0.677395 
0.677395 
0.819439 
0.819439 


321 168 429.685 
18 061876 000 
303 016 377.265 20 
23 914 008 800 


0.010 
150 000 
0.000 11 
150 000 



From theoretical considerations and from experimen- 
tal results with Rb, a quadratic i?-dependence of A is 
expected. If we assume that A{B) = A^) x [1 + kB"^), 
then there arc five unknowns to solve for: Aq, k, gi / g.j. 
Xi, and X2. In addition to the four equations for the 
low magnetic field measurements, derived from the Breit- 
Rabi formula, a fifth equation is given by the expression 
for the high-B value of A: 

Ao X (1 + kBs^) = -625 008 837.371 Hz, (9) 

where Bs = 4.4609 T. Solving the set of five equations 
gives 

Ao = -625 008 837.044(12) Hz, (10a) 
k = 2.63(18) X 10"" (10b) 
gj'/gj = 2.134 779 852 0(23) x lO"''. (10c) 

The uncertainties of the parameters are obtained by vary- 
ing the experimental frequencies through their uncertain- 
ties. The value of gi' / gj given by Eq. (10c) is consistent 
with, but less precise than, that obtained from the high- 
B data alone [Eq. (6)]. 



iih. electron. The interaction with an external mag- 
netic field B — Bz is taken into account by the mini- 
mal coupling prescription, i.e., making the replacement 
Pi ^ Pi + eA(ri), where A is the vector potential func- 
tion, Vx A = B. Nuclear Zeeman, electron spin Zeeman, 
and hyperfine interactions are also added as perturba- 
tions to Hq . The kinetic energy term for the ith electron 
in Eq. (11) undergoes the change 



PL 

2m 



2A(r,)]^ 



2m 

Pf. , f [Pi • A(ri) + A(ri) • pj 

2m 

+ 



2m 



2m 



EL 

2m 



Hf + Ht 



(12a) 

(12b) 
(12c) 



The term containing the first power of A is called the 
paramagnetic interaction Hj, while the one containing 
A^ is called the diamagnetic interaction Hf. The divi- 
sion into paramagnetic and diamagnetic parts is gauge- 
dependent, but the choice of gauge, 



IV. CALCULATION OF DIAMAGNETIC 
HYPERFINE SHIFT COEFFICIENT k 

In nonrelativistic atomic theory, the diamagnetic shift 
in hyperfine structure arises as a cross term involving 
both the diamagnetic interaction and the hyperfine in- 
teraction in second-order perturbation theory. Let the 
unperturbed Hamiltonian for an A/'-electron atom with 
nuclear charge Ze be 



N 2 
1=1 



N 

y — 

^ 47reorj 
1=1 



E 



47reo|ri 



(11) 



A(r) = X B, 



(13) 



is particularly convenient. With that choice, the param- 
agnetic term becomes 



/fP = — [p,-A(r,)+A(r,)-p,] 

= —{Pi • X B -h Tj X B • Pj 
4to 



(14a) 

(14b) 



e ^ e ^ eUi)zB ,^ , , 

= 2;;^^^^p^-^ = 2;;^^-^ = ^' ^^^^^ 

whore li = r; x p.^ is the orbital angular momentum of 
the ith electron. The diamagnetic term becomes 



where m, is the electron mass, — e is the electron charge, 
and Ti and Pj are the position and momentum of the 



2m 



8m 



(rixB)^ 



8to 



ix'i+yf)B^. (15) 
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Hf can be divided into a spherically symmetric (scalar) 
part and a rank- 2 spherical tensor part: 

„2o2 2 „2o2 2 

Hf = £^^ + £,^(l-3cos2 0O = Hf+Hf. (16) 



12m 



24m 



The total paramagnetic interaction is obtained by sum- 
ming Hf over all electrons: 



i=l 



2m 



(17) 



where is the z-component of the total electronic or- 
bital angular momentum. To a very good approxima- 
tion, the ground electronic state is an S'-state, that is, an 
eigenstate of with L = 0, so ff p can be neglected. 
The total scalar part of the diamagnetic interaction is 



N 
i=l 



N 2d 

X - e B 



i=l 



2„2 



12m 



(18) 



In first-order perturbation theory, this leads to a com- 
mon shift of all of the hyperfine-Zeeman sublevels of the 
groimd electronic state. In second-order perturbation 
theory, there is a cross term that is first-order in both 
iJ'^" and the Fermi contact hyperfine interaction. This 
leads to an efi^ective interaction that appears as a shift in 
the hyperfine A value, proportional to B"^ [7, 8, 29-31]. 
Evaluation of the second-order perturbation term yields 
the constant k. 

Unpublished calculations by Lipson, similar to those 
done for Rb [9, 31], yielded k = 2.52 x 10"" T-^ for 
Be"*". These calculations used Hartree-Fock-Slater wave 
functions [32]. The inclusion of all intermediate states, 
including the continuum, was done by solving an inho- 
mogeneous differential equation for the perturbed wave 
function [33, 34]. Similar calculations were done by one 
of the present authors (W.M.I.) , but using a paramet- 
ric potential for Be+ that reproduces the experimental 
energy levels [35]. This yielded k = 2.68 x 10"" T'^. 
The estimate given in Ref. [5] of /S.A = -0.017 Hz T-^ 
(equivalent to fc = 2.7 x 10~^^ T~^) was based on these 
calculations. 

Another method of obtaining k to the same order in 
perturbation theory is to calculate an approximate elec- 
tronic wave function that is acciiratc to first order in 
H^'^ and then to calculate the mean value of the hy- 
perfine interaction with these wave functions. A sim- 
ple way to do this is to use the MCHF (multiconfig- 
uration Hartree-Fock) or MCDHF (multiconfiguration 
Dirac-Hartree-Fock) method, where an additional term, 
6r? is added to each single-electron Hamiltonian. Un- 
like model potential methods, MCHF and MCDHF are 
ah initio in the sense that they require no experimen- 
tal input, such as observed energy levels, only values of 
fundamental constants. 

The GRASP set of MCDHF programs [36-39] were 
used to calculate correlated wave functions for Be^ with 



and without the hrf term and to calculate hyperfine con- 
stants with these wave functions [40]. The calculation 
of the unperturbed hyperfine constant for ^Be"*" was 
similar to that done by Bieron et al., [41] but less exten- 
sive, leaving out for example nuclear recoil and the Breit 
interaction. The result, = —624.19 MHz, is within 
0.13 % of the experimental value. The calculation was 
then modified by including the hrf terms in the Hamil- 
tonian. In Hartree atomic units (e = m = 7i = 1), 6 is 
dimensionless. It was varied from 1 x 10~^ to 2 x 10~^. 
The change in A, relative to Aq, was found to be propor- 
tional to h for 5 < 1 X 10~^. Evaluating the constant of 
proportionahty yields k — 2.645(2) x 10^^^ T^^, where 
the uncertainty here reflects only numerical error, not er- 
ror due to physical approximations, such as neglect of the 
Breit interaction. However, given the good agreement of 
the calculated and experimental values of Aq, we estimate 
the error of the calculated value of k to be no more than 
1 %. All of the calculated values of k, including Hartree- 
Fock-Slater, parametric potential, and MCDHF, arc in 
good agreement with the experimental result within the 
experimental error of 7 %. 

The use of the diamagnetic potential in a relativistic 
calculation requires some justification, since the minimal 
coupling prescription does not yield a term proportional 
to in the Dirac Hamiltonian [42]. Instead of Eq. (12a) 
we have, for the kinetic energy term in the single-electron 
Dirac Hamiltonian: 

ca • Pj ca • [pj -I- eA(rj)] = ca • Pj -|- eca • A(ri), (19) 

where 



a 



cr 
cr 



(20) 



and (T = ((T""^, cr^, (7^) is defined in terms of the Pauli ma- 
trices cr*. Since the field-dependent perturbation only 
contains B to the first power, calculation of k with 
Eq. (19) requires third-order perturbation theory (two 
orders in the magnetic field interaction and one order in 
the hyperfine interaction) , unlike the nonrelativistic case, 
which requires only second-order perturbation theory. 

Kutzelnigg [43] showed that a unitary transformation 
of the Dirac Hamiltonian in the presence of a magnetic 
field yields terms resembling the nonrelativistic paramag- 
netic and diamagnetic terms, plus another term which is 
first-order in A and whose effect goes to zero in the non- 
relativistic limit. This justifies the use of second-order 
perturbation theory to calculate k in the relativistic case. 
The relativistic form of the single-electron diamagnetic 
interaction is [42, 43] 



(21) 



which differs from the nonrelativistic form [Eq. (15)] only 
by the factor of /3, where ^ is the 4x4 matrix 



/ 
-/ 



(22) 
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where / is a 2 x 2 identity matrix. The factor of /3 
was found also by Szmytkowski by a different method 
[44]. The effect of /3 on a matrix element of Hf{rcl) is 
to reverse the sign of the integral involving the product 
of the small components of the Dirac orbitals. Thus, 
the relative error inciirred by ignoring 8 should be less 
than {Za)^, where a here is the fine-structure constant 
e^/(47reo?ic). This error can be neglected for Be+ {Z = 4) 
but may amount to a few percent for Rb (Z = 37). 

The total tensor part of the nonrelativistic diamagnetic 
interaction is 



JV 



N 



H 



d2 



E^f = E^(l-3cos^^.)- (23) 

i=l j=l 



In first-order perturbation theory, this leads to no energy 
shifts in the ground electronic state. In second-order per- 
turbation theory, there is a cross term that is first-order 
in both H'^^ and the electric quadrupolc hyperfine in- 
teraction. This leads to an effective interaction called 
the magnetically induced quadrupole hyperfine interac- 
tion [9, 31]. Comparison of the second-order perturba- 
tion expression for the indiiccd quadrupole interaction 
and the expression for the quadrupole antishielding fac- 
tor 7oo defined by Sternheimer [45] yields 



fQ 



2Amh 



(24) 



where /q is the coefficient of the induced quadrupole 
interaction defined in Eq. (7). This form is useful be- 
cause values of have already been calculated for 
many atoms and ions. The value 700 = 0.7088 for 



Be"*" has been calculated in a Hartree-Fock approxima- 
tion [46]. The most recent experimental value of the 
^Be nuclear quadrupole moment is 5.288(38) x 10"^" m^ 
[47]. These values of the constants yield the estimate 
/q = 6.64 X 10~® Hz T~^, which is much smaller than 
the experimental upper limit set in Sec. HE. 

The hyperfine-assisted Zeeman shift can be estimated 
by the same method as that used for the rubidium atom 
[11]. In this approximation, the hyperfine matrix ele- 
ments are given by the Fermi- Segre formula [48], the ns- 
state energies are obtained from a hydrogenic (quantum 
defect) approximation, and the continuum s-state wave 
functions are obtained from a Coulomb approximation. 
In this approximation, the coefficient /?hz for ^Be+ is 
equal to 2.61 x 10~^ Hz T~^, which is smaller than the 
experimental upper limit set in Sec. HE. 

Some other B-dependent shifts are in principle present, 
such as a magnetic-field-induced spin-dipole hyperfine 
term, but based on calculations done for Rb [31] they are 
likely to be much smaller than the terms already consid- 
ered. Various corrections to gj and to g/ (e.g., nuclear 
diamagnetic shielding) are calculable but are beyond the 
scope of this paper. For recent calculations of gj for Be"*" 
and other three-electron atoms, see Refs. [49-52]. For a 
recent calculation of the nuclear diamagnetic shielding of 
Be+, see Ref. [53]. 
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